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Abstract 

We discuss leading charm production in connection with energy deposition in 
the central rapidity region. Special attention is given to the correlation between 
production in central and fragmentation regions. If the fraction of the reaction 
energy released in the central region increases the asymmetry in the xp distribu- 
tions of charmed mesons will become smaller. We illustrate this quantitatively with 
simple calculations performed using the Interacting Gluon Model. Leading beauty 
production is also considered. 
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Several experiments have reported |jl| a significant difference between tlie xp dependence 
of leading and nonleading cliarmed mesons. It was not possible to explain these data with the 
usual perturbative QCD [|| or with the string fragmentation model contained in PYTHIA and 
some alternative mechanisms have been advanced. The most detailed data analysis, including 
predictions for asymmetry and leading particle effect at higher energies, has been carried out with 
the intrinsic charm model (ICM). In this model the essential ingredient for a good description 
of asymmetries is the recombination mechanism which binds together the intrinsic (fast) charm 
quarks and the valence quarks of the projectile. Apart from this fast component there is a 
slow one, which populates predominantly the central rapidity (low xp) region and is given by 
perturbative QCD. The ICM is thus a two-component model where the central (parton fusion) 
and fragmentation regions (containing intrinsic charm) components are completely independent 
and added in a simple way. In particular, there is no energy conservation constraint imposed on 
both components, which would obviously result in some simple kinematical correlations between 
them. 

The purpose of this work is to show that such kinematical correlation between central and 
non-central production is relevant for the study of the observed asymmetries in the produc- 
tion of charmed mesons and that it is also connected to another characteristic of high energy 
multiparticle production processes, namely to the inelasticity K of the reaction. It defines the 
fraction of the initial energy y/s which is released and deposited in the central region of reaction. 
In particular its energy (-y/i) dependence will be important here. All models that adress charm 
production in the central region predict that there is no asymmetry in this region. Asymmetry 
comes from the fragmentation (large rapidity y) region. Therefore, if K increases with energy 
there will be less energy available in the large y region and this will result in a softening of the 
leading xp distributions. Notice that this is independent of all ingredients of the hadronization 
process since they are universal and energy independent (like, for example, the fragmentation 
functions). The xp distributions of the leading charmed particles will thus eventually merge with 
the distribution of the centrally produced charmed particles, which will then become broader. 
The asymmetry will then not be observed! The opposite might also be true if K decreases with 
energy. In this case the leading system will carry proportionally more and more energy, implying 
a faster leading charm and resulting in stronger asymmetry 
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The asymmetry problem can therefore be weh formulated just in terms of kinematical consid- 
erations. All dynamics will show itself only in the way through which initial energy of projectiles 
will be distributed in rapidity space. For example, one would naively expect that, if perturbative 
QCD becomes more important at higher energies (due to, for example, increased minijet activ- 
ity), the central production (and also energy deposition in central region) will become dominant 
and the asymmetry will decrease or even disappear. This goes along with the expected increase 
with energy of inelasticity deduced from the analysis of accelerator and cosmic ray data @. 

In what follows we shall therefore study leading charm production in terms of the Interacting 
Gluon Model (IGM) ||^, which has been invented to describe the inelasticity and its energy 
behaviour and recently used also to successfuly describe many aspects of multiparticle production 
(including its semi-hard minijets component, which can be important for charm production at 
high energies). 

Since the IGM has already been described previously in great detail Q we shall provide 
here, for completeness, only the most basic formulas and concentrate our attention on the 
specific mechanisms of charm production and on the calculation of the asymmetries between 
leading and nonleading charm mesons. The asymmetry has been most accurately measured in 
the np scattering, therefore we shall start discussing this process first Q. In Fig. la we show the 
IGM description of the energy flow in a hadron- hadron collision at high energies. Through the 
cooperative action of a certain number of soft gluons, carrying an overall momentum fraction xi 
of the incoming pion, colliding with a similar bunch of gluons coming from the target nucleon 
and carrying fraction X2 of its momentum, an object called central fireball (CF) is formed. It 
will decay later on producing observed secondaries. In the IGM the probability for this to 
happen is given by the function xixi,X2)- The pion remnants leaving the central region (i.e., 
their valence quarks plus some gluons which did not interact) carrying momentum xl are called 
in the IGM the leading jet (LJ) and, being themselves excited objects, may also produce particles 
(including D mesons). 
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From the basic function xixi,X2) we can compute the Feynman momentum distributions of 
the CF, x{xcf), where xcf = xi — X2, and of the LJ, fuixi), by a simple change of variables: 



where m£, (1.8 GeV) and mo are the masses of the D meson and of the lightest state produced 



fragmentation production. The momentum distributions of the systems which will later give 
origin to charmed particles are derived from the same quantity x{xi-,X2)- Moreover, x{xcf) 
and fLjixi) carry all the energy (\/s) dependence of the process, which is both explicit, in the 
theta functions, and implicit, since x{xi,X2) depends on ^/s . In Fig. lb we show central DD 
meson production where D{D) is any D meson carrying a c(c) quark. Notice that, in the spirit 
of IGM, the central production ignores the valence quarks of target and projectile (defined here 
as those which carry the essential quantum numbers of the colliding pions and protons) which, 
in the first approximation, just "fly through". Because of this, the centrally produced D's will 
not show any leading particle effect. There are two distinct ways to produce D mesons out 
of LJ's: fragmentation and recombination. It is assumed here that, whenever energy allows, 
we shall have also cc pairs in the LJ (produced, for example, from the remnant gluons present 
there). These charmed quarks may undergo fragmentation into D mesons, as shown in Fig. Ic, 
but may also as well recombine with the valence quarks as depicted in Fig. Id. It turns out 
that only this last process will produce asymmetry. In the case of pion-nucleon scattering, the 
measured leading charmed mesons are D~ and the nonleading are D'^. 




(1) 





in such collisions In the above equations we clearly see the connection between central and 
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We shall now write the Feynman xp single inclusive distribution of D mesons produced by 
the CF, by the fragmentation in the LJ (F) and by the recombination there (R): 

^ — = / dxcF x{xcf) I dxc g{xc) D f-^) , (3) 



dx 



I dxL fLj{xL) I dxc g{xc) D f-^) , (4) 



D- •'^D- •'^n- ^ Xc 

da^ fi 



dx 



dXL fLj{xL) dxc dxc dxu dxd g{xc) g{xc) f{Xu) f{Xd) 
D- ■^^D- J J J J 

■ 6{xj:)- - Xc- Xd)S{xL - Xc- Xc- Xd- Xu) , (5) 

where f{x) and g{x) are distribution functions of valence and charm quarks respectively and the 
D{z)^s are charm quark fragmentation functions The momentum distributions are given 
by @ and (Q) with the replacements: D^, c ^ c. These nonleading mesons will not be 

produced by recombination (eq. (|5|)). The functions f{x) and g{x) are essentially unknown since 
they are momentum distributions of partons inside the CF and inside the LJ after the collision. 
Following our assumption that the valence quarks interact weakly we shall approximate f{x) by 
the initial state valence quark distributions and take them from reference Q. As for the charm 
quark distribution, g{x), the situation is less clear. The c — c pairs do not come directly from the 
sea : in the CF they are produced and in the LJ they may be excited. It is therefore reasonable 
to think that the charm quarks will be somewhat faster than ordinary sea quarks. Accordingly 
we shall use for g{x) the ansatz proposed by Barger and collaborators [p^]: 

. ^ fi-xy/^ 

g{x) = (6) 



which is less singular than 1/x but still much softer than an intrinsic charm distribution which 



behaves typically like x{\ — x). The fragmentation functions have the Peterson form |11 



N 

Dc^d{z) = -7- — 

z[l — 1/z — — z)f 

where e ~ -, — , and m„, pqT, mo, Pqt are mass and transverse momentum of the light and 
of the heavy quark respectively and is a normalization constant. In the present case £ ~ 0.06. 
In Fig. 2a we show the (unnormalized) contributions coming from the three processes above (eqs. 
@ and (|5|)). As expected, central production (solid line) leads to the softest D meson xp 



distribution, recombination in the leading jet (dotted line) leads to the hardest final distribution 
and leading jet fragmentation (dashed line) lies in between. This is so because xi^cp) is softer 
than fLji^L} and because recombination adds momenta whereas fragmentation causes always 
some deceleration. Note that, although flat, the dashed and dotted curves have a pronounced 
maximum at very low xq. This is a direct consequence of the behaviour of g{x). If instead of 
the form (6) we use an intrinsic charm distribution we will obtain a strong suppression at low x 
and a maximum around xo = 0.4 — 0.6 . A final comment on this figure is that our distribution 
of centrally produced D's (solid line) is broader than the one obtained from perturbative QCD. 
This is so because the cooperative mechanism adds together soft gluons, increasing the energy 
released in the central region, favouring higher values of xi and X2 (in Fig. la) and allowing for 
fluctuations with higher xcf- Considering all that was said above we can conclude that the IGM 
(like the ICM) is a two component model in which the components are not very much different 
in shape from each other (in sharp contrast to what happens with the components of the ICM) 
and have some overlap. Because of this we expect to find smaller asymmetries than those found 
in ref. but this depends, of course, on how one mixes the different components. In what 
follows we write the differential cross section as the sum of a central fireball (CF) and leading 
jet (LJ) component and the last one as the sum of a fragmentation (F) and a recombination (R) 
component, using a similar notation as in ref. Q: 

a ~ d^ + ^ d^ ' 

1 da^-^ , 1 da^ 1 da^ 

(i-O-px— + (9) 



a^-^ dx£)~ dxf)- dx^, 

where the mixture parameters are ^ (0 < < 1) and 

^ = ,CF + ,LJ (10) 

In the case of the D"^ distribution, the expressions above are the same but ^ = 0. 

The ICM parameter rj has been chosen 0.2 because of an analogy between aic (our cr^"^) 
and the diffractive charm cross section. On the other hand, in the Valon Model the same 
data are adressed without any central component. This would correspond to taking r/ = 1. 
Here, because of the kinematical mixing between CF and LJ the value of r/ is essentially free. 
In what follows we will choose it to be ry = 0.7. Note also that, in our case, ^ = corresponds 
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to no asymmetry. Since existing data on open charm production [T| apparently do not show 
nuclear effects |13], we use here (as all other models which address these data) the IGM for 
hadron-nucleon collisions [14|. 



In Fig. 3 we compare our calculations with WA82 data. Fig. 3a and 3b show the xp 
spectrum of and , respectively, and Fig. 3c shows the asymmetry which is given by: 

dtr da 

Ai^F) = ^ (11) 



+ 



dx 



D+ 



In Fig. 3b and 3c solid, dashed and dotted lines correspond to ^ = 0.8, 0.5 and 0.2 respec- 
tively. Data points are from the WA82, E769 and E791 collaborations. As it can be seen, 
a satisfactory description of data can be obtained with the IGM. The best description can be 
obtained with a large ammount of recombination (^ = 0.8). This is ultimately due to our choice 
of g{x). We have checked that the choice of an ordinary sea distribution for the charm quarks 
in the CF and LJ requires r/ = 1.0 and = 1.0 for a reasonable fit. The choice of an intrinsic 
charm distribution allows for small values of rj and ^. The conclusion seems to be that although 
data do not rule out usual sea distributions as an input, good fits with more reasonable values 
of the parameters can be obtained using harder charm quark distributions like (6) or the one 
used in ref. 4. 

We consider now the energy dependence of the asymmetry. All details concerning the par- 
ticularities of charm production are energy independent. In equations (8) and (9) r], ^ and the 
differential distributions, i.e., respectively normalization and shape of the curves, can depend on 
■y/s. For simplicity we shall assume that ^ does not change with the energy. The distributions 
will depend on ^/s through xi^cp) and fLj{xL)- The behaviour of these last functions 
with \/s is shown in Fig. 4a and 4b respectively. We observe a very modest broadening of 
xi^cp) implying a small increase of (xcp) and a more pronounced softening of /lj(2^l) with 
the corresponding reduction of (xl)- As for rj, an extensive analysis of charged particle pro- 
duction up to Tevatron energies has shown that it decreases by a factor 3 when we go from ISR 
to Tevatron energies. Assuming a similar reduction for the case of charmed particle production 
T] will change from 0.7 to 0.25. Considering what was said above we evaluate again all the 
expressions (l)-(9) at y/s = 1800 GeV. The resulting asymmetry is shown in Fig. 4c with a 
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dashed line. For comparison we show in the same figure (with a sohd hne) the asymmetry at 
^/s = 26 GeV calculated with the same parameters. It decreases 20% in the region xp > 0.5. 
Although this is not a very impressive change it illustrates the trend. Moreover, we know that 
the asymmetry goes asymptotically to zero since 77 ^ 0. We emphasize that this is so because 
the IGM (in its version ||^]) predicts that at higher energies, because of the action of minijets, 
the energy deposition in the central region will increase implying two effects: a growth of the 
central multiplicity ( implying thus an increase of 1 — r/ ) and a softening of the leading jet 
momentum distribution. We can therefore conclude that, irrespective of details of charm pro- 
duction, these both effects combined will reduce the asymmetry. It is interesting to mention 
that the data collected in Fig. 3c come from three different collaborations E769, WA82 and 
E791 with beam energies of 250, 340 and 500 GeV respectively. In the CMS this corresponds 
to a variation of y/s = 23 GeV to ^/s = 33 GeV. This energy change is small, the error bars 
are large and therefore no change in the asymmetry is visible yet. At higher energies there is a 
chance to experimentally verify this behaviour at RHIC or LHC. 

As a straightforward extension of our analysis we calculate now the asymmetry in B meson 
production. This is done by simply replacing m£, by niB = 4.75 GeV in (1) and e = 0.06 by 
£ = 0.006 in (7). In principle we should also change g{x) but in a first estimate we keep the 
ansatz (6). If we would use an intrinsic distribution for g{x) it would be very weakly dependent 
on the heavy quark mass Q. In Fig. 5 we show the xi^cp) distribution for charm (dashed 
line) and for beauty production (solid line) with the proper change in eq. (1). The energy 
is ^/s = 26 GeV. The effect of increasing the production threshold {m£, ^ in the theta 
function in eq.(l)) is to select events with a more massive CF and with larger lower limits for 
xi and X2i suppressing thus larger values of xci;' = xi — X2 with respect to charm production 
(in the limit of total energy deposition, i.e., xi = X2 = 1, the CF would be at rest). This 
effect is however very small. This is expected and seen in Fig. 5. In Fig. 6a (6b) we show 
the Xp distributions of nonleading (leading)Z) and B mesons. The energy is the same as in 
Fig. 5 and the parameters are the same as before (r/ = 0.7 and ^ = 0.8). Nonleading spectra 
are calculated with eqs. (3) and (4). The Peterson fragmentation functions, appearing in those 
equations, are very sensitive to the value of e. In the case of beauty, the strong reduction of 
e makes the fragmentation function strongly peaked at very large values of z. The emerging 
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B^s will therefore be much less decelerated than the -D's. This effect compensates the previous 
one and the final nonleading B distribution is harder than the nonleading D one. The leading 
distributions include recombination, given by eq. (5), which is not affected by the change in the 
heavy quark mass. Because of this, the spectra in Fig. 6b exhibit the same qualitative behaviour 
(i?'s faster than D's) seen in Fig. 6a but the difference between i?'s and D^s is smaller. The 
asymmetries of B^ /B^ and / are shown in Fig. 7 with solid and dashed lines respectively 
for ^/s = 26 and 1800 GeV. The asymmetry in beauty is about 50% weaker than in charm at 
xp = 0.8 and both show a similar decrease with energy. 

In conclusion: the IGM describes the energy flow in high energy hadron collisions. It takes 
properly into account the correlation between energy deposition in the central region and the 
leading particle momentum distribution. It accounts for charmed meson production in a natural 
and satisfactory way and makes the prediction that at higher energies the increase of inelasticity 
K (see 0]) will lead to the decrease of the asymmetry in heavy quark production. It also predicts 
a weaker asymmetry for beauty. We believe that this point should also be adressed by other 



models which deal with asymmetry in heavy flavour production [16|. 
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Figure Captions 

Fig. 1 (a) Illustration of a pion-nucleon collision: fractions xi and X2 of the incoming hadrons 
momenta form a central fireball (CF) with probability x(xi, X2). The fraction 1 — xi = xl 
is carried by the leading jet (LJ). The leading jet momentum spectrum is fuixL)', (b) 
Nonleading D meson production by central fireball fragmentation; (c) Nonleading D meson 
production by leading jet fragmentation; (d) Leading D meson production by leading jet 
recombination . 

Fig. 2 xp distributions of D mesons calculated with the IGM. The solid line shows the contri- 
bution of the central fireball fragmentation. The dashed line shows the contribution of 
the leading jet fragmentation and the dotted line shows the contribution of the leading 
jet recombination. 

Fig. 3 (a) xp distribution of mesons calculated with our model and compared with WA82 
data; (6) xp distribution of D" mesons calculated with our model and conpared with 
WA82 data. Solid line corresponds to ^ = 0.8 in eq. (9) while dashed and dotted lines 
correspond to ^ = 0.5 and ^ = 0.2 respectively; (c) the asymmetry calculated with the 
IGM and compared with WA82 (solid circles), with E769 (open squares) and E791 (open 
triangles) data. Solid, dashed and dotted lines correspond to the same choices of ^ of (3b). 

Fig. 4 (a) Momentum distribution xi^cp) of the central fireball at ^/s = 26 GeV (solid line) 
and at s/s = 1800 GeV (dashed line) ; (6) momentum distribution fuixL) of the leading 
jet at ^/s = 26 GeV (solid line) and at s/s = 1800 GeV (dashed line); (c) the asymmetry 
in pion- proton collision calculated with the IGM: solid line corresponds to ^/s = 26 GeV 
with T] = 0.7 and dashed line corresponds to ^/s = 1800 GeV with 7] = 0.25. In both cases 
e = 0.8. 

Fig. 5 xp distribution {x{xcf)) of centrally produced b — b (solid line) and c — c (dashed line) 
quark pairs. 

Fig. 6 (a) xp distribution of nonleading B (solid line) and D (dashed line) mesons; (6) the same 
as (a) for leading B (solid line) and D (dashed line) mesons. The energy is in both cases 
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^/s = 26 GeV. 

Fig. 7 B~ /B^ (solid lines) and D~ /D^ (dashed lines) asymmetries at ^/s = 26 and 1800 GeV. 
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